Comprehensive proteogenomic characterization of rare kidney tumors.

Non-clear cell renal cell carcinomas (non-ccRCCs) encompass diverse malignant and benign tumors. Refinement of differential diagnosis biomarkers, markers for early prognosis of aggressive disease, and therapeutic targets to complement immunotherapy are current clinical needs. Multi-omics analyses of 48 non-ccRCCs compared with 103 ccRCCs reveal proteogenomic, phosphorylation, glycosylation, and metabolic aberrations in RCC subtypes. RCCs with high genome instability display overexpression of IGF2BP3 and PYCR1. Integration of single-cell and bulk transcriptome data predicts diverse cell-of-origin and clarifies RCC subtype-specific proteogenomic signatures. Expression of biomarkers MAPRE3, ADGRF5, and GPNMB differentiates renal oncocytoma from chromophobe RCC, and PIGR and SOSTDC1 distinguish papillary RCC from MTSCC. This study expands our knowledge of proteogenomic signatures, biomarkers, and potential therapeutic targets in non-ccRCC.

Systemic Amyloidosis and Kidney Transplantation: An Update.

Amyloidosis is a heterogeneous disorder characterized by abnormal protein aggregate deposition that often leads to kidney involvement and end-stage kidney disease. With advancements in diagnostic techniques and treatment options, the prevalence of patients with amyloidosis requiring chronic dialysis has increased. Kidney transplantation is a promising avenue for extending survival and enhancing quality of life in these patients. However, the complex and heterogeneous nature of amyloidosis presents challenges in determining optimal referral timing for transplantation and managing post-transplantation course. This review focuses on recent developments and outcomes of kidney transplantation for amyloidosis-related end-stage kidney disease. This review also aims to guide clinical decision-making and improve management of patients with amyloidosis-associated kidney disease, offering insights into optimizing patient selection and post-transplant care for favorable outcomes.

Proteogenomic data and resources for pan-cancer analysis.

The National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC) investigates tumors from a proteogenomic perspective, creating rich multi-omics datasets connecting genomic aberrations to cancer phenotypes. To facilitate pan-cancer investigations, we have generated harmonized genomic, transcriptomic, proteomic, and clinical data for >1000 tumors in 10 cohorts to create a cohesive and powerful dataset for scientific discovery. We outline efforts by the CPTAC pan-cancer working group in data harmonization, data dissemination, and computational resources for aiding biological discoveries. We also discuss challenges for multi-omics data integration and analysis, specifically the unique challenges of working with both nucleotide sequencing and mass spectrometry proteomics data.

Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation.

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Single cell characterization of myeloma and its precursor conditions reveals transcriptional signatures of early tumorigenesis.

Multiple myeloma is a plasma cell malignancy almost always preceded by precursor conditions, but low tumor burden of these early stages has hindered the study of their molecular programs through bulk sequencing technologies. Here, we generate and analyze single cell RNA-sequencing of plasma cells from 26 patients at varying disease stages and 9 healthy donors. In silico dissection and comparison of normal and transformed plasma cells from the same bone marrow biopsy enables discovery of patient-specific transcriptional changes. Using Non-Negative Matrix Factorization, we discover 15 gene expression signatures which represent transcriptional modules relevant to myeloma biology, and identify a signature that is uniformly lost in abnormal cells across disease stages. Finally, we demonstrate that tumors contain heterogeneous subpopulations expressing distinct transcriptional patterns. Our findings characterize transcriptomic alterations present at the earliest stages of myeloma, providing insight into the molecular underpinnings of disease initiation.

Immune biomarkers of response to immunotherapy in patients with high-risk smoldering myeloma.

Patients with smoldering multiple myeloma (SMM) are observed until progression, but early treatment may improve outcomes. We conducted a phase II trial of elotuzumab, lenalidomide, and dexamethasone (EloLenDex) in patients with high-risk SMM and performed single-cell RNA and T cell receptor (TCR) sequencing on 149 bone marrow (BM) and peripheral blood (PB) samples from patients and healthy donors (HDs). We find that early treatment with EloLenDex is safe and effective and provide a comprehensive characterization of alterations in immune cell composition and TCR repertoire diversity in patients. We show that the similarity of a patient's immune cell composition to that of HDs may have prognostic relevance at diagnosis and after treatment and that the abundance of granzyme K (GZMK)+ CD8+ effector memory T (TEM) cells may be associated with treatment response. Last, we uncover similarities between immune alterations observed in the BM and PB, suggesting that PB-based immune profiling may have diagnostic and prognostic utility.

Host protein kinases required for SARS-CoV-2 nucleocapsid phosphorylation and viral replication.

Multiple coronaviruses have emerged independently in the past 20 years that cause lethal human diseases. Although vaccine development targeting these viruses has been accelerated substantially, there remain patients requiring treatment who cannot be vaccinated or who experience breakthrough infections. Understanding the common host factors necessary for the life cycles of coronaviruses may reveal conserved therapeutic targets. Here, we used the known substrate specificities of mammalian protein kinases to deconvolute the sequence of phosphorylation events mediated by three host protein kinase families (SRPK, GSK-3, and CK1) that coordinately phosphorylate a cluster of serine and threonine residues in the viral N protein, which is required for viral replication. We also showed that loss or inhibition of SRPK1/2, which we propose initiates the N protein phosphorylation cascade, compromised the viral replication cycle. Because these phosphorylation sites are highly conserved across coronaviruses, inhibitors of these protein kinases not only may have therapeutic potential against COVID-19 but also may be broadly useful against coronavirus-mediated diseases.

High-Resolution Profiling of Lung Adenocarcinoma Identifies Expression Subtypes with Specific Biomarkers and Clinically Relevant Vulnerabilities.

Lung adenocarcinoma (LUAD) is one of the most common cancer types and has various treatment options. Better biomarkers to predict therapeutic response are needed to guide choice of treatment modality and to improve precision medicine. Here, we used a consensus hierarchical clustering approach on 509 LUAD cases from The Cancer Genome Atlas to identify five robust LUAD expression subtypes. Genomic and proteomic data from patient samples and cell lines was then integrated to help define biomarkers of response to targeted therapies and immunotherapies. This approach defined subtypes with unique proteogenomic and dependency profiles. Subtype 4 (S4)-associated cell lines exhibited specific vulnerability to loss of CDK6 and CDK6-cyclin D3 complex gene (CCND3). Subtype 3 (S3) was characterized by dependency on CDK4, immune-related expression patterns, and altered MET signaling. Experimental validation showed that S3-associated cell lines responded to MET inhibitors, leading to increased expression of programmed death-ligand 1 (PD-L1). In an independent real-world patient dataset, patients with S3 tumors were enriched with responders to immune checkpoint blockade. Genomic features in S3 and S4 were further identified as biomarkers for enabling clinical diagnosis of these subtypes. Overall, our consensus hierarchical clustering approach identified robust tumor expression subtypes, and our subsequent integrative analysis of genomics, proteomics, and CRISPR screening data revealed subtype-specific biology and vulnerabilities. These LUAD expression subtypes and their biomarkers could help identify patients likely to respond to CDK4/6, MET, or PD-L1 inhibitors, potentially improving patient outcome. SIGNIFICANCE: Integrative analysis of multiomic and drug dependency data uncovers robust lung adenocarcinoma expression subtypes with unique therapeutic vulnerabilities and subtype-specific biomarkers of response.

Molecular map of chronic lymphocytic leukemia and its impact on outcome.

Recent advances in cancer characterization have consistently revealed marked heterogeneity, impeding the completion of integrated molecular and clinical maps for each malignancy. Here, we focus on chronic lymphocytic leukemia (CLL), a B cell neoplasm with variable natural history that is conventionally categorized into two subtypes distinguished by extent of somatic mutations in the heavy-chain variable region of immunoglobulin genes (IGHV). To build the 'CLL map,' we integrated genomic, transcriptomic and epigenomic data from 1,148 patients. We identified 202 candidate genetic drivers of CLL (109 new) and refined the characterization of IGHV subtypes, which revealed distinct genomic landscapes and leukemogenic trajectories. Discovery of new gene expression subtypes further subcategorized this neoplasm and proved to be independent prognostic factors. Clinical outcomes were associated with a combination of genetic, epigenetic and gene expression features, further advancing our prognostic paradigm. Overall, this work reveals fresh insights into CLL oncogenesis and prognostication.

Genetic subtypes of smoldering multiple myeloma are associated with distinct pathogenic phenotypes and clinical outcomes.

Smoldering multiple myeloma (SMM) is a precursor condition of multiple myeloma (MM) with significant heterogeneity in disease progression. Existing clinical models of progression risk do not fully capture this heterogeneity. Here we integrate 42 genetic alterations from 214 SMM patients using unsupervised binary matrix factorization (BMF) clustering and identify six distinct genetic subtypes. These subtypes are differentially associated with established MM-related RNA signatures, oncogenic and immune transcriptional profiles, and evolving clinical biomarkers. Three genetic subtypes are associated with increased risk of progression to active MM in both the primary and validation cohorts, indicating they can be used to better predict high and low-risk patients within the currently used clinical risk stratification models.

Single-nucleus cross-tissue molecular reference maps toward understanding disease gene function.

Understanding gene function and regulation in homeostasis and disease requires knowledge of the cellular and tissue contexts in which genes are expressed. Here, we applied four single-nucleus RNA sequencing methods to eight diverse, archived, frozen tissue types from 16 donors and 25 samples, generating a cross-tissue atlas of 209,126 nuclei profiles, which we integrated across tissues, donors, and laboratory methods with a conditional variational autoencoder. Using the resulting cross-tissue atlas, we highlight shared and tissue-specific features of tissue-resident cell populations; identify cell types that might contribute to neuromuscular, metabolic, and immune components of monogenic diseases and the biological processes involved in their pathology; and determine cell types and gene modules that might underlie disease mechanisms for complex traits analyzed by genome-wide association studies.

ZBTB33 is mutated in clonal hematopoiesis and myelodysplastic syndromes and impacts RNA splicing.

Clonal hematopoiesis results from somatic mutations in cancer driver genes in hematopoietic stem cells. We sought to identify novel drivers of clonal expansion using an unbiased analysis of sequencing data from 84,683 persons and identified common mutations in the 5-methylcytosine reader, ZBTB33, as well as in YLPM1, SRCAP, and ZNF318. We also identified these mutations at low frequency in myelodysplastic syndrome patients. Zbtb33 edited mouse hematopoietic stem and progenitor cells exhibited a competitive advantage in vivo and increased genome-wide intron retention. ZBTB33 mutations potentially link DNA methylation and RNA splicing, the two most commonly mutated pathways in clonal hematopoiesis and MDS.

A proteogenomic portrait of lung squamous cell carcinoma.

Lung squamous cell carcinoma (LSCC) remains a leading cause of cancer death with few therapeutic options. We characterized the proteogenomic landscape of LSCC, providing a deeper exposition of LSCC biology with potential therapeutic implications. We identify NSD3 as an alternative driver in FGFR1-amplified tumors and low-p63 tumors overexpressing the therapeutic target survivin. SOX2 is considered undruggable, but our analyses provide rationale for exploring chromatin modifiers such as LSD1 and EZH2 to target SOX2-overexpressing tumors. Our data support complex regulation of metabolic pathways by crosstalk between post-translational modifications including ubiquitylation. Numerous immune-related proteogenomic observations suggest directions for further investigation. Proteogenomic dissection of CDKN2A mutations argue for more nuanced assessment of RB1 protein expression and phosphorylation before declaring CDK4/6 inhibition unsuccessful. Finally, triangulation between LSCC, LUAD, and HNSCC identified both unique and common therapeutic vulnerabilities. These observations and proteogenomics data resources may guide research into the biology and treatment of LSCC.

The FDA-approved drug Alectinib compromises SARS-CoV-2 nucleocapsid phosphorylation and inhibits viral infection in vitro.

While vaccines are vital for preventing COVID-19 infections, it is critical to develop new therapies to treat patients who become infected. Pharmacological targeting of a host factor required for viral replication can suppress viral spread with a low probability of viral mutation leading to resistance. In particular, host kinases are highly druggable targets and a number of conserved coronavirus proteins, notably the nucleoprotein (N), require phosphorylation for full functionality. In order to understand how targeting kinases could be used to compromise viral replication, we used a combination of phosphoproteomics and bioinformatics as well as genetic and pharmacological kinase inhibition to define the enzymes important for SARS-CoV-2 N protein phosphorylation and viral replication. From these data, we propose a model whereby SRPK1/2 initiates phosphorylation of the N protein, which primes for further phosphorylation by GSK-3a/b and CK1 to achieve extensive phosphorylation of the N protein SR-rich domain. Importantly, we were able to leverage our data to identify an FDA-approved kinase inhibitor, Alectinib, that suppresses N phosphorylation by SRPK1/2 and limits SARS-CoV-2 replication. Together, these data suggest that repurposing or developing novel host-kinase directed therapies may be an efficacious strategy to prevent or treat COVID-19 and other coronavirus-mediated diseases.

Genomic Profiling of Smoldering Multiple Myeloma Identifies Patients at a High Risk of Disease Progression.

PURPOSE: Smoldering multiple myeloma (SMM) is a precursor condition of multiple myeloma (MM) with a 10% annual risk of progression. Various prognostic models exist for risk stratification; however, those are based on solely clinical metrics. The discovery of genomic alterations that underlie disease progression to MM could improve current risk models. METHODS: We used next-generation sequencing to study 214 patients with SMM. We performed whole-exome sequencing on 166 tumors, including 5 with serial samples, and deep targeted sequencing on 48 tumors. RESULTS: We observed that most of the genetic alterations necessary for progression have already been acquired by the diagnosis of SMM. Particularly, we found that alterations of the mitogen-activated protein kinase pathway (KRAS and NRAS single nucleotide variants [SNVs]), the DNA repair pathway (deletion 17p, TP53, and ATM SNVs), and MYC (translocations or copy number variations) were all independent risk factors of progression after accounting for clinical risk staging. We validated these findings in an external SMM cohort by showing that patients who have any of these three features have a higher risk of progressing to MM. Moreover, APOBEC associated mutations were enriched in patients who progressed and were associated with a shorter time to progression in our cohort. CONCLUSION: SMM is a genetically mature entity whereby most driver genetic alterations have already occurred, which suggests the existence of a right-skewed model of genetic evolution from monoclonal gammopathy of undetermined significance to MM. We identified and externally validated genomic predictors of progression that could distinguish patients at high risk of progression to MM and, thus, improve on the precision of current clinical models.

Scaling computational genomics to millions of individuals with GPUs.

Current genomics methods are designed to handle tens to thousands of samples but will need to scale to millions to match the pace of data and hypothesis generation in biomedical science. Here, we show that high efficiency at low cost can be achieved by leveraging general-purpose libraries for computing using graphics processing units (GPUs), such as PyTorch and TensorFlow. We demonstrate > 200-fold decreases in runtime and ~ 5-10-fold reductions in cost relative to CPUs. We anticipate that the accessibility of these libraries will lead to a widespread adoption of GPUs in computational genomics.

Prediction of pKa Using Machine Learning Methods with Rooted Topological Torsion Fingerprints: Application to Aliphatic Amines.

The acid-base dissociation constant, pKa, is a key parameter to define the ionization state of a compound and directly affects its biopharmaceutical profile. In this study, we developed a novel approach for pKa prediction using rooted topological torsion fingerprints in combination with five machine learning (ML) methods: random forest, partial least squares, extreme gradient boosting, lasso regression, and support vector regression. With a large and diverse set of 14 499 experimental pKa values, pKa models were developed for aliphatic amines. The models demonstrated consistently good prediction statistics and were able to generate accurate prospective predictions as validated with an external test set of 726 pKa values (RMSE 0.45, MAE 0.33, and R2 0.84 by the top model). The factors that may affect prediction accuracy and model applicability were carefully assessed. The results demonstrated that rooted topological torsion fingerprints coupled with ML methods provide a promising approach for developing accurate pKa prediction models.

SNPs2ChIP: Latent Factors of ChIP-seq to infer functions of non-coding SNPs.

Genetic variations of the human genome are linked to many disease phenotypes. While whole-genome sequencing and genome-wide association studies (GWAS) have uncovered a number of genotype-phenotype associations, their functional interpretation remains challenging given most single nucleotide polymorphisms (SNPs) fall into the non-coding region of the genome. Advances in chromatin immunoprecipitation sequencing (ChIP-seq) have made large-scale repositories of epigenetic data available, allowing investigation of coordinated mechanisms of epigenetic markers and transcriptional regulation and their influence on biological function. To address this, we propose SNPs2ChIP, a method to infer biological functions of non-coding variants through unsupervised statistical learning methods applied to publicly-available epigenetic datasets. We systematically characterized latent factors by applying singular value decomposition to ChIP-seq tracks of lymphoblastoid cell lines, and annotated the biological function of each latent factor using the genomic region enrichment analysis tool. Using these annotated latent factors as reference, we developed SNPs2ChIP, a pipeline that takes genomic region(s) as an input, identifies the relevant latent factors with quantitative scores, and returns them along with their inferred functions. As a case study, we focused on systemic lupus erythematosus and demonstrated our method's ability to infer relevant biological function. We systematically applied SNPs2ChIP on publicly available datasets, including known GWAS associations from the GWAS catalogue and ChIP-seq peaks from a previously published study. Our approach to leverage latent patterns across genome-wide epigenetic datasets to infer the biological function will advance understanding of the genetics of human diseases by accelerating the interpretation of non-coding genomes.